Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN I*TeX style file v2.2) 



PSR B1828— 11: a precession pulsar torqued by a quark 
planet? 



K. Uu\ Y. L. Yue^, and R. X. Xu^ 

^Department of Geophysics, School of Earth and Space Science, Peking University, China 
^Department of Astronomy, School of Physics, Peking University, China 



o 
o 

(N 



5 February 2008 



(N 
> 

(N 



o 
o 



ABSTRACT 

The pulsar PSR B1828— 11 has long-term, highly periodic and correlated variations 
in both pulse shape and the rate of slow-down. This phenomenon may provide evi- 
dence for precession of the pulsar as suggested previously within the framework of free 
precession as well as forced one. On a presumption of forced precession, we propose 
a quark planet model to this precession phenomenon instead, in which the pulsar is 
torqued by a quark planet. We construct this model by constraining mass of the pulsar 
(Mpsr), mass of the planet (Mpi) and orbital radius of the planet (rpi). Five aspects 
are considered: derived relation between Mpsr and Tpi, movement of the pulsar around 
the center of mass, ratio of Mpsr and Mp\, gravitational wave radiation timescale of the 
planetary system, and death-line criterion. We also calculate the range of precession 
period derivative and gravitational wave strength (at earth) permitted by the model. 
Under reasonable parameters, the observed phenomenon can be understood by a pul- 
sar (10~^ ~ 10~^Mq) with a quark planet (10~® ^ IO^^Mq) orbiting it. According to 
the calculations presented, the pulsar would be a quark star because of its low mass, 
which might eject a lump of quark matter (to become a planet around) during its 
bhth. 

Key v^rords: Pulsars: individual (PSR B1828— 11) — stars: planetary systems — 
gravitational waves 



1 INTRODUCTION 

The pulsar PSR B 1828— 11 shows long-term, highly pe- 
riodic and correlated variations in both the pulse shape 
and the slow-down rate. Its variations are best described 
as harmonically related sinusoi ds, with periods of approxi- 
mately 1000, 500 and 250 days l|Stairs et a l.' 2000). The phe- 
nom enon indicates the mos t compelling evidence for preces- 
sion l|Link fc Epsteinll200ll ). 

To explain this phenomenon, some authors 
(IJones fc AnderssonI 1200 ll : [r ezanial |2003| ) have proposed 
different models within the framework of free precession. 
The observation could not be a problem in the standard 
vie w of neutron st a rs if t he star's crust is free to precess. 
In iLink fc EpsteinI (|200ll ). the correlated changes in the 
pulse duration and spin period derivative can be explained 
as a precession of the star's rigid crust c oupled to the 
magnetic dipole torque. lAkgun et al.l (|2006l ) modelled the 
timing behavior with the inclusion of both geometrical 
and spin-down contributions to the residuals. However, 
investigations concerned on internal structure of neutron 
stars show that free precession may be damped out if 
vortices pinning to the stellar crust and hydrodynamics 
(MHD) coupling are taken into consideration. In detail, the 



rotation of the superfluid, accounting for a large proportion 
of the moment of inertia of the pulsar, is contained in 
an array of vortices. Models, in which vortices pinned to 
the stellar crust become unpinned during a glitch, might 
have described the occurrence of and recovery from glitches 
(|Alpar et al.l ll984l V The vortex pinning will damp out 
free prec ession on time s cales of several hundre d precession 
periods l|Shahaml 1 19771 : ISedrakian et al.1 Il999l ) if the pin- 
ning force is as strong as suggested in the glitch models. 
Additionally, the MHD coupling between the crust and 
the core will also strong ly affect precession of the pulsar 
(I Levin fc D 'AngelollioM ). The decay of precession, caused 
by the mutual friction between the neutron superfluid and 
the plasma in the core, is expected to occur over tens to 
hundreds of precession periods an d may be measu rable 
over a human lifetime. As noted by iLinkI l|2003l . [200^ . the 
picture of vortex lines entangled in flux tubes appears to be 
incompatible with observations of long-period precession, 
which indicates that the standard scenario of the outer 
core (superfluid neutrons in co-existence with type II, 
superconducting protons) should be reconsidered. 

An alt ernative w ay is to consider the pulsar as a solid 
quark star l|Xull2003l ). where precession models will not re- 
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ceive puzzle that damping out brings. But there are still 
some problems when here we come to the model of free pre- 
cession. For example, the ellipticity (or dynamical flatten- 
ing) of the pulsar derived from free precession model is not 
fitted well with the one calculated by Maclaurin approx- 
imation. Consider the pulsar as a rotation ellipsoid with 



the principal moment of inertia 



ly < Iz and the 



corresponding radii a — b > c. In free precession mod- 
els, the stellar dynamical flattening is e = (/^ — Ix)/Ix ~ 
e^/(2 - e^) = P/Pproco ~ 10"**, where P is the spin 
period, Pprocc is the precession period of the pulsar and 
e = A^/l — c? jo? is the stellar eccentricity. It can be approx- 
imated as since e ^ 1. Meanwhile, from Maclaurin 
approximation, the stellar ellipticity can be determined by 
e = \l-(c/af [']/(c/af^^ ^eV s « 26/3 « 3x lO'^ProL « 
2 X 10"^ (|Xull2006l : [Zhou et all [20041 '). These two values, e 
and e, which are expected to be generally matched if free 
precession model is available, are quite different. Therefore, 
new ideas need to be devised to explain the phenomenon 
of precession instead of the free precession ones. Actually, 
a forced precession model driven by an fossil disk was pre- 
sented in lQiao et all (|2003t )F1 

Here we alternatively present a quark planet model to 
explain the phenomenon of precession. In this model, forced 
precession is caused by a quark planet orbiting the pulsar. 
In Sect. 2, flrst we establish the relation between mass of 
the pulsar and orbital radius of the planet in case that the 
dynami cal flatten ing is obtained from Maclaurin approxi- 
mation l|Xull2006l ). Then we explain why the planet should 
be a quark planet, rather than a normal one like the earth 
or Jupiter when planet model is referred to. Next we limit 
movement of the pulsar around the center of mass by er- 
rors in the TO As (time-of-arrival). Death- line criterion and 
limitation on gravitation wave radiation timescale are also 
considered so as to constrain orbital radius, mass of the pul- 
sar and mass of the planet. In Sect. 3, we calculate the pre- 
cession period derivative and gravitational wave radiation 
strength of the pulsar for different mass of the pulsar and 
orbital radius of the planet. In Sect. 4, we conclude by dis- 
cussing the formation of such system and expecting further 
observation to test the model. 



2 PRECESSION TORQUED BY A QUARK 
PLANET 

First of all, we suppose that the pulsar PSR1828— 11 could 
be either neutron star or quark star since both are candi- 
date models for pulsar. In case that the precession period 
is much more than the spin period and the orbital period, 
the angular ve locity of forced precession can be expressed as 
iMenke fc Ab bott 2004,) 



3GA/pi 

— r— ecos t^pi. 



(1) 



^ In this paper, the authors obtained the stellar oblateness from 
Maclaurin approximation and considered the precession as the 
whole star's motion, not only the crust's. So we think that here 
an idea of solid quark star is better than neutron star so as to 
prevent decay of precession. 
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Figure 1. Relation between mass of the planet (Afpi) and or- 
bital radius (rpi) for 500-day precession period from Eq. (2). We 
can have Mpi cos 0pi fa Mpi if 9pi is not close to 90° . The figure 
indicates that reasonable value of Mp (Mpi < Mq) may be found 
while Tpi is less than 10^ cm. 



where a is the precession angular velocity, u) is the rotational 
angular velocity of the spinning pulsar, is orbital radius 
of the planet, G is gravitational constant, Mpi is mass of 
the planet, e ~ 3e/2 ~ 3 x 10^'' is the stellar dynamical 
flattening and 6'pi is average inclination of planet orbit. From 
Eq. (1) we can have 



8^^ 



3GPPp 



-''pi, 



(2) 



where P = 27r/a'p is spin period and Pprecc = 2n /a is pre- 
cession period of the pulsar. If we consider a normal planet 
similar to the earth or Jupiter, the typical value of rpi should 
be 0.1 AU or 1 AU and the corresponding value of A/pi is 
much larger than the solar mass. In Fig. 1, the relation be- 
tween Mp cos Spi and Tpi is shown derived from the 500-day 
precession period. We can see that if rpi reaches 10^ cm and 
9pi is not close to 90° , mass of the planet will be over the 
solar mass. A planet needs to be of several billion times of 
Mq to provide enough torque if it locates at 1 AU away from 
the pulsar. We cannot believe the existence of such a planet 
since it would definitely induce huge orbital timing effects in 
the pulse residuals. However, the result is not surprising be- 
cause the pulsar has a much shorter forced precession period 
and thus the torque that dominates the precession needs to 
be much stronger. Meanwhile, precession torque is reduced 
when the distance between the pulsar and the planet be- 
comes longer. Consequently, if there is a planet close to the 
pulsar, it may be able to provide enough torque to cause the 
short-period precession. That is the reason that we consider 
quark planet since its orbital radius could mainly depend on 
the kick energy, which can vary in a large range (Sect. 4). 

Therefore, we suppose that the orbital radius of the 
planet is between (10^ cm, 10^ cm), where 10*' cm is typical 
radius of a normal neutron star. Besides, for PSR B1828— 11, 
the errors in the TOAs (time-of-a rrival) are limited by ran- 
dom noise to about Tc ~ 0.2 ms (|Stairs et al.ll2000l) . So in 
such a planetary system, the pulsar is not likely to move 
more than about TcC « 6 x 10^ cm around the center of mass, 
and its orbital radius should be less than 3 x 10^ cm. In ad- 
dition, if the eccentricity of the orbit is not considerable, we 
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Figure 2. Constraint on rpi (tlie orbital radius of the planet) 
and Mpsr (the mass of the pulsar) by observations and theoretical 
arguments. The shadowed 'Available' region surrounded by Line 
(1)— (5) is the parameter space for rpi and Mpsr. The five lines are 
defined by Eq. (20)— (24). In this figure we use cosSpi = 1, since 
the locations of Line (l)-(3) will not change much with variation 
of average inclination of planet orbit 8pi from 0° to 80° (see Table 



can have MpiVpi « Mpsrrpsr, where Mpsr and rpsr are mass of 
the pulsar and its orbital radius around the center of mass, 



respectively. Since rpsr has a maximum rp 
the relation can be derived as Mpirpi < M, 



,=3 X 10" cm, 
Fi- 



psr ' psr,max 



nally, in such a planetary system, mass of the pulsar should 
be much larger than that of the planet so we approximately 
assume Mpar/Mpi > fc = 10. 

Next we consider the gravitational wave radiation 
(GWR) of the planetary system for further limitation. In 
normal double neutron star system, the distance between 
the two stars is about 10^'' cm. Ther efore the timescale of 
GWR is rather long usua lly 10* years (|Hulse fc Tavlorll 19751 : 
iTavlor fc WeisbereJl 19821 ). However, in this quark planet sys- 
tem, due to the short distance between the two objects, the 
power of GWR may be significantly larg er. In double-star 
system it is given by (|Misner et al.|[l973 ) 



dE _ 320" /i'm^ _ 32G*MgsrMpVMpsr + MpQ 
dt ^ 



(3) 



5c^ 5c^a^ 
where a ~ rpi is the semi-major axis of the orbit, m = Afpi -|- 
Mpsr and /i = MpiMpsr/(Afpi -|- Mpsr) is the reduced mass. 
The GWR costs the total of potential energy and dynamic 
energy of the planetary system 

GMplMpsr 



-Etot = — - 



2r 



pi 



So timescale of GWR can be derived as 

IStotI 



(4) 



(5) 



dE/dt 64G3 MpiMpsr (Mpsr + Mpi) • 

If the planetary system is stable, the timescale must be long 
enough. Here we approximately set it as r > ro = 10* 
yearsRi 3 x 10^^ s. 

Additionally, we consider the death-line criterion, which 
requests the potential drop at the polar cap of the pulsar 
be more than (j)o ~ 1 0^^ V l|Ruderman fc Sutherland!! 19751 : 
!Usov fc Melrosel!l995! '). If we assume that PSR1828-11 is an 
aligned pulsar, potential drop is 



ttWB . 2 
:r — sm ( 

cP 



(6) 



where B is the polar magnetic field strength at pulsar 
surface, R ~ (3Mpsr/(47rp))^''^ is the pulsar radius, p ~ 
7 X 10^* g/cm'^ is the density of the pulsar and 6 = 
arcsin ^J2-kR/{cP) is the opening half-angle of the polar cap. 
Here we use the density for quark stars to obtain the lower 
limit of Mpsr. The magneti c field can be approximated by 
!Manchester fc Tavloj 



B 



Mc^PP 



(7) 



where I ~ (2/5)Mpsri?^ is the principal moment of inertia. 
From Eq. (6) and (7), the relation between potential drop 
and mass of the pulsar can be derived as below 



,3,i/2,37r^,l/3, P sl/2/ 1n1/3 ,.,5/6 



(8) 



While > 00, we have 



20000^^^^^ 



-^y/^ « 3 X lo-^M©. 



(9) 



Mpsr > ( 

2437r4p3 

Actually, the assumption of alignment in PSR1828— 11 is 
rather strong. The potential drop from Eq. (6) can be larger 
by more than one order of magnitude if the inclination of th e 
magnetic axis to the spin axis is not zero (!Yue et al.!!2006! ). 
Consequently, constraint on mass of the pulsar can be lower 
by about one magnitude and thus we have Ms > 10"*Mq. 
Now there are five limitations for Mpsr, '"pi and Mpi: 



Mpsr/Mpl > k, 

Mpirpi < Mpsrr psr, max ; 



5cv:^ 



> To, 



64G3MpiMpsr(Ms-hMpi) 

rpi € (10® cm, lO'' cm). 



Mpsr > 10"* M( 



0- 



(10) 
(11) 

(12) 

(13) 
(14) 



If we consider Afpsr 2> Mpi and substitute for A4pi in term of 
rpi according to Eq. (2), then the limitations can be derived 
as 

S'^-' ^ (15) 



Mp, 



-'"pi, 



M 



> 



8^^ 



psr ' 3PPprcccGeC0Sepirp; 

15PP, 



4 

-'"pi, 



sr .max 



Mpsr < ( 



prccc"- 



-''pi) 



1/2 



5127r2G2rocos Opi 
10® cm < rpi < 10^ cm, 



Mpsr > IO'^Mg 



(16) 

(17) 

(18) 
(19) 



In Fig. 2 we consider the above limitations and figure 
out the available range for Mpsr and rpi. Accordingly, point 
(rpi, Mpsr) should be above Line (1) and (2), below Line 
(3) and (5) and on the right of Line (4). So we have the 
shadowed area as the "Available" area for point (rpi, Mpsr). 
Line (l)-(5) are defined as below 



Line (1) : Mp- 
Line (2) : Mp, 



8fc7r'' 



8^^ 



3PPp rece Ge cos 6pirp, 



sr. max 



-'"pi, 



4 

-^pi. 



(20) 
(21) 
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Line (3) : Mp^r = ( 



15PPn 



5127r2G2rocos( 



Line (4) : rpi 
Line (5) : Mp^r ~ 



10" cm, 



0- 



From Fig. 2, the available value range of Vpi, 

, (10- 

9piM0 



and Mpi are (10** cm, 5 x lO'^ cos^/^ Sp 
6 X lO-^cos-^''^6lpiM0) and (6 x 10" 



cmj 



(22) 

(23) 
(24) 

Mpsr 

, 6 X 



10" cos" ' 6piMq). Here we use cos6pi = 1 because the 
positions of Line (1), (2) and (3) do not vary distinctly with 
the changing of 9pi from 0° to 80° . Different value ranges of 
Tpi, Mpar and Mpi with different 9pi are shown in Table 1. 
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3 TO TEST THE MODEL BY FURTHER 
OBSERVATION 

Loss of the total energy of the system caused by gravita- 
tional wave radiation will lead to decay of the planet or- 
bit. Correspondingly, precession period will be reduced with 
the decreasing of orbital radius. Meanwhile, the planetary 
system may act as a detectable gravitational wave source. 
Therefore, it is possible to test and improve the model by 
GWR detection and long-period observation for precession 
period derivative (Pproco). 

So next we calculate Pprccc and the characteristic am- 
plitude of GWR source strength (/ic) for different rpi and 
Mpar, respectively. From Eq. (2), if orbital radius has a slight 
change, the variation of Pproco can be expressed as 



APp; 



GPAfpiecosSpi^^'P'- ^^^^ 

In this case we do not consider the change of spin period 
(The result will prove its reasonableness). Similarly, from 
Eq. (4) we have 

GMpiMp, 



^»^pl 



-Arpi, 



(26) 



In addition, the slight change of mechanical energy is caused 
by GWR in a short period (see Eq. [3]) 

,p 32G^Mp\rMp^(Mpa. + Mpi) ^^ 

" '^d^^ ^ ' 

Considering Mpar 3> Mp\ and combining Eq. (25)-(27) give 
the period derivative of precession as below 



APp, 



At 



5c^ Pecos QpxTpi 



(28) 



Meanwh ile, rate of loss of angu lar momentum caused by 
GWR is lUshomirskv et all l|2000h 



(29) 



where E is the rate of loss of the total energy, 57 — 
^/GMpar/r^i is the period of revolution of the planet and 
d=3.58 kpc (|Tavlor fc CordesI [l993l ') is the distance of the 
pulsar, /la is the source's 'angle- averaged' field strength (at 
earth) and approx i matel y we have ha ~ he {ha ~ 1.15/ic, see 
lUshomirskv et al.1 (|200G| ')'). Combining Eq. (2), (3) and (29) 
gives 



h 



3V5dc4PPprece e cos 6pi 



M, 



psr ' pi- 



(30) 



Figure 3. A zoomed parameter space for the 'Available' re- 
gion. The procession period derivative (Pprece = dPprece/dt, dash 
lines), the perturbed metric (ha, dot lines), and gravitational wave 
frequency (u, dash-dot lines) are drawn. Here we use cos 9pi = 1 
in the calculations. The available area gives model-permitted pa- 
rameter space for Pprccc, hs, and u. 



Besides, the frequency of GWR is 



27r 



GMp 



3 

'^pi 



(31) 



In Fig. 3, relations between r-pi and Mpi from Eq. (28) 
for a group of Pprccc, from Eq. (31) for a group of h^ and 
from Eq. (32) for a group of u are shown (cos6lpi ~ 1). The 
relations are limited by the available area for point (rpi, Mpi) 
from Fig. 1. As is shown, the maximum and minimum of 
Pprccc are 4x lO""* and 6 x 10"'' while those of /la are 3x 10~^^ 
and 1 X 10"^" . The result indicates that the precession period 
changes much quickly than spin period of the pulsar and 
GWR at earth is not intense enough to be detected by LIGO 
at its working frequency. For example, at a frequency of 10 
Hz, value of /la is about lO"'^^, which is below the current 
detection limit of the LIGO at the same frequency (about 
10-22). 



4 CONCLUSION AND DISCUSSION 

Within the framework of forced precession, we propose 
a quark planet model to explain the precession of PSR 
B1828— 11. The observed phenomenon can be understood 
by a pulsar (probably a quark star) together with a quark 
planet which torques dominantly the pulsar to precess. In 
principle, orbital radius of the quark planet should be be- 
tween lO'' cm and 10** cm while the range of mass of the pul- 
sar and the planet are approximately (lO-'^M©, 10~^Mq) 
and (IQ-^M©, 10"^ Mq), respectively. These resuhs might 
not be strange since other candidate s of low-m ass q uark 
stars were also discussed previously (|Xul l2005l : lYue et al.l 
l2006h . We calculate the model-permitted precession period 
derivative and characteristic amplitude of GWR for the sys- 
tem. The precession period changes much quickly than spin 
period of the pulsar; meanwhile, GWR strength at earth 
may not be large enough to be detected by current LIGO. 

If there is a quark planet providing torque for the forced 
precession of pulsar PSR B 1828— 11, it should be close to the 
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Table 1. The parametric range of rpi, Afpsr and Mpi for different inclination of planet orbit, dpi. The 
variation of the range is not significant for 9pi from 0° to 80° . 



6»pi Tp Ms Mp 



0° 


(lO^* 


cm, 5 X lO'^ 


cm) 


(10- 


"*M0, 6 X IQ- 




(6 X 10" 


-^Mq, 6 X 10" 




30° 


(lO^* 


cm, 5 X 10^ 


cm) 


(10- 


-*Mq, 6 X IQ- 


-3M0) 


(7 X 10- 


-^Mq, 7 X 10" 


-*Mq) 


60° 


(106 


cm, 5 X 10^ 


cm) 


(10- 


-*Mq, 8 X IQ- 


-''Me) 


(1 X 10- 


■*Mq, 9 X 10" 


-*Mq) 


80° 


(106 


cm, 4 X 10'' 


cm) 


(10- 


"^Mq, 1 X IQ- 


-2M0) 


(3 X 10- 


■^Mq, 2 X 10" 


-SMq) 



pulsar with a distance of several times of the pulsar's radius. 
The pulsar mass should also be significantly lower than Mq , 
which may suggest that the pulsar would be a quark star. 
Such kind of planets, orbiting closely to the center pulsars, 
could be ejecta during the formation of the quark stars with 
stron g turbulence if the surface energy is reasonably low (|Xul 
I2OO6I ). Considering the orientation of the system's angular 
momentum, the planet is not likely to have a inclination of 
orbit very close to 90° and our previous analysis with 6'pi 
varying from 0° to 90° can work effectively. 

In this paper we do not consider the possibility of more 
than one planet, which may provide a way to explain the 
other two possible precession periods of the pulsar. The lim- 
itation on orbital radius and ratio of the pulsar mass to the 
planet one could be improved if the formation of the sys- 
tem is considered. However, precession periods of the pulsar 
cannot be exactly obtained now from the seemingly periodic 
post-fit timing residuals. Long-period timing observations in 
the future are necessary in order to obtain more accurate 
precession period derivative of the system. Since 2000, the 
pulsar has accomplished several precession period. There- 
fore, if precession period derivative reaches its maximum in 
this model, the precession period may have changed several 
days. In addition, we need observation for gravitational wave 
to test and improve the model. Whether it can be detected 
or not will both provide further limitation on mass of the 
pulsar and orbital radius of the planet. 

The model could also be tested by X-ray observation. 
(1) If the pulsar is a solid quark star with Mpar ~ 1Q~'Mq 
and then rpsr ~2 km, the rate of rotation energy loss could 
be only about 10^ times smaller than that in the standard 
model where the pulsar is a normal neutron star. Assum- 
ing that only ~ 0.1% of the spin-dow n power could turn 
into the non-thermal X-ray luminosity (jLorimer fc Kramer! 
I2OO5I ). the fiux at earth should be about 1 x 10 erg ■ 
cm"^ • s"^ (less than 1 photon/60 hours). (2) The ther- 
mal X-ray emissivity from the pulsar could also be lower. 
If thermal emission is from the global star, the flux should 
be ~ 5x 10"^^erg-cm-^-s"^ and ~ 3x 10"" erg-cm-^-s"^ 
for the pulsar with surface temperatures of 200 eV and 100 
eV, respectively. Taking absorption into consideration, we 
can expect a fiux of ~ lO"^* erg • cm"^ • s"^ (about 70 pho- 
tons/6 hours). But if the pulsar is a neutron star with radius 
10 km and surface temperature > 60 eV0 the flux is much 
higher, > 10"^^ erg • cm"^ • s"^. Future observations of the 
pulsar by Chandra or XMM-Newton could certainly bring 
us more details about the real nature. 



According to the standard cooling model by IPag j l|l998l ). the 
temperature of this ~ 10^ years old pulsar is ^ 68 eV. 



Finally, we note that the nature of pulsars (to be neu- 
tron or quark stars) is still a matter of debate even after 
40 years of the discovery. The reason for this situation is in 
both theory (the uncertainty of non-perturbative nature of 
strong interaction) and observation (the difficulty to distin- 
guish them). It is a non- mainstream idea that pulsars are 
actually quark stars, but this possibility cannot be ruled 
out yet according to either first principles or observations. 
"Low-mass" is a natural and direct consequence if pulsars 
are quark stars since quark stars with mass < IMq are 
self-confined by color interaction rather by gravity. An ar- 
gument against the low-mass idea could be the statistical 
mass-distribution of pulsars in binaries (~ lA Mp,). How - 
ever this objection might not be so strong due to l|Xull2005l ): 
(1) if the kick energy is approximately the same, only solar- 
mass pulsars can survive in binaries as low-mass pulsars may 
be ejected by the kick; (2) low-mass bare strange stars might 
be uncovered by re-processing the timing data of radio pul- 
sars if the pulsars' mass is not conventionally supposed to 
be ~ I.4M0 . In this work, we just try to understand the pe- 
culiar precession nature of PSR B1828-11 in the quark star 
scenario, since the mainstream-scientific solution to preces- 
sion might not be simple and natural. 
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